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We have performed a theoretical study of the magnetic circular dichroism in the x-ray ab- 
sorption spectra (XMCD) of the equiatomic CoPt and FePt ordered alloys as a function of the 
spin quantization axis. We found that the magnetization axis is along the [001] direction and the 
magneto-crystalline anisotropy energy (MCA) for the FePt compound is twice as large as that of the 
CoPt compound in agreement with experiment. The band structure and the total density of states 
confirm that all electronic states contribute to the MCA, and not just the states at the vicinity of 
the Fermi level. The orbital magnetic moments decrease with respect to the angle between the [001] 
axis and the spin quantization axis, and are much larger for the CoPt compound. We show that the 
00 , orbital moment anisotropy is reflected in the XMCD signal. 
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I. INTRODUCTION 



^ . The disordered equiatomic binary alloys,of the XY (X= Fe, Co - Y= Pd, Ft) type crystallize in the fee strueture and 
the magnetization is along the [111] axis.EJ At low temperatures these alloys tend to order in the LIq layered-ordered 
structure and in this case the spontaneous magnetization tends to align perpendicular to the layer stacking explaining 
the behavior of CoFt films during the magnetic annealingfl The strong perpendicular magnetic anisotropy (FMA) is 
due to the highly anisotropic LIq structure and it makes them very attractive for magnetic recording devices.El The 
Q ■ latter structure can be also obtained by molecular beam epitaxy (MBE) of alternating layers of pure X and Y atoms 
Q due to the substrate induced constraints. The first observation of the Llg long-range order for a [001] CoFt film grown 
by MBE was made by Harp et aUa in 1993 and for a [001] FeFt film by CeboUada et alB Lately other techniques have 
been also employed to develop films presenting FMA. CoFt films were grown by sputtering by Viisokay et ala and by 
evaporation by Lin and Gorman,td and FeFt films were grown by various sputtering techniques .Bllj 

■ The magneto-crystalline anisotropy energy (MCA) can be probed by many techniques such as torque or ferro- 
T— I \ magnetic resonance measurements. Both these methods-|describe the MCA in terms of phenomenological anisotropy 

■ constants. It has been demonstrated by Weller et aLll3 that x-ray magnetic circular dichroism (XMCD) is also a 
suitable technique for probing the MCA, via the determination of the anisotropy of the orbital magnetic moment on 
a specific shell and site. The x-ray absorption spectroscopy (XAS) using polarized radiatkuipprobes element specific 
magnetic properties of alloys by applying the XMCD sum rules to the experimental spectra.llJ"E3 However for itinerant 



systems, in pa:|^g^ to low symmetry systems, the use of the sum rules is debated because they are derived from an 



' atomic theory.tjOllj Lately angle-dependent XMCD experiments Jiave been used to provide a deeper understanding 
] for the relation between MCA and the orbital magnetic moments.E3 ■— ■ 
I • Thejc-ray absorptipa for CoFt multi-layejrs has been alreadv studied experimentally by Nakajima et al.f^ Koide 
' et al.f^ R-uegg et al.fn and Schiiltz et a/E3 Makajima et a/E3 revealed a strong enhancement of the cobalt orbital 
Ch moment when FMA was present, Koide et alVA showed that with decreasing cobalt thickness the easy axis rotates 
^ from in-plane to outj^-plane, and Riiegg et alzA that platinum polarization increases also with decreasing cobalt 
thickness. Hlil et alc^ showed by x-ray absorption spectroscopy that modifications of platinum edges in different 
compounds are correlated to the change in the number of holes. ___ 

Several ab-initio calculations have already been performed to investigate the XMCD.E3'cSE3 The L2- and Ls-edges 
involving electronic excitations of 2p-core electrons towards d- valence states have primarily attracted much attention 
d I due to dependence of the dichroic spectra on the exchange-splitting and the spin-orbit coupling of both initial core 
and final valence states. For 5d elements dissolved in 3d transition metals, the spin-orbit coupling of the initial 2p-core 
states is large and the resulting magnetic moment is small, while the opposite is true for the Zd elepijents. This can 
lead to a pronounced dichroic spectra as seen by Schiitz in the case of 5d elements dissolved in iron.c2l 

In this work we study the correlation between the quantization axis dependent XMCD and the magnetic properties 
of both ordered alloys FeFt and CoFt. Ouc_rpethod is based on an all-electron relativistic and spin-polarized full- 
potential muffin-tin orbital method (LMTO^E3E3 in conjunction with both the von Barth and HedipjDarameterization 
to the local density approximation (LSDA)Eil and the generalized gradient approximation (GGA)E3 to the exchange 
correlation potential. The implementation of the calculation of the XMCD spectra has been presented in a previous 
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work.c3 In section 2 we present the details of the calculations, while in sections 3 and 4 we discuss our MCA and the 
magnetic spin and orbital moments, respectively. In section 5 we present our calculated XMCD as a function of the 
spin quantization axis, and in section 6 we discuss the interpretation of the MCA using the band structure and the 
total density of states anisotropy. 

II. COMPUTATIONAL DETAILS 

To compute the ejlfictronic properties of CoPt and FePt jsze used the experimental lattice constants (a=3.806A and 
c/a=0.968 for CoPtO; a= 3.86lA and c/a=0.981 for FePta) and a unit cell containing one atom of cobalt(iron) and 
one of platinum. The LIq structure can be seen as a system of alternating cobalt (iron) and platinum layers along the 
[001] direction. The MCA and the XMCD are computed with respect to the angle 7 between the [001] axis and the 
spin quantization axis on the (010) plane. So 7=0° corresponds to the [001] axis and 7=90° to the [100] axis. 

MCA can be computed directly using ab-initio methods; it is defined as the difference between the total energy 
for two different spin quantization axis. The spin-orbit coupling contribution to MCA is implicitly included in om 
ab-initio calculations, and we do not take into acpaunt the many-body interactions of the spin-magnetic momentsE^ 
since their contribution to the MCA is negligibleEj The number of k points for performing the Brillouin zone (BZ) 
integration depends strongly on the interplay between theucontributions to the MCA from the Fermi surface and the 
remaining band structure contribution to the total energy.Ej When the former contribution to the MCA is important, 
a large number of k-points is needed to describe accurately the Fermi surface. For the two studied systems we found 
that 6750 k-points in the BZ are enough to converge the MCA within 0.1 mcV. 

To perform the integrals over the BZ we use a Gaussian broadening method which convolutes each discrete eigenvalue 
with a Gaussian function of width 0.1 eV. This method is known to lead to a fast and stable convergence of the spin 
and charge densities compared to the standard tetrahedron method. To develop the potential inside the MT spheres 
we calculated a basis set of lattice harmonics including functions up to ^ = 8, while for the FFT we used a real space 
grid of 16x16x20. We used a double set of basis functions, one set to describe the valence states and one for the 
unoccupied states. For the valence electrons we used a basis set containing 3xs, 3xp and 2xc? wave functions, and 
for the unoccupied states 2xs, 2xp and 2xd wave functions. 

III. MAGNETO-CRYSTALLINE ANISOTROPY 

CoPt and FePt films are known to present a strong uniaxial MCA, because of the highj^isotropic LIq inter-metallic 
phase. Experimentally the magnetization axis is found to be along the [001] direction.ESO In a first step we performed 
calculations with 250 k-points in the BZ. This number of k-points is large enough to produce accurate total energy 
when Gaussian smearing is used for the integration in the Brillouin Zone but not enough accurate to compute the 
MCA. In figure |l| we present calculations with respect to the angle 7 between the [001] axis and the spin quantization 
axis on the (010) plane for the CoPt compound within the LSDA. We observe that total energy value increases with 
the angle 7. So the ground state corresponds to 7=0°, i.e. the [001] axis. The same behavior occurs for the FePt 
compound. Because 250 k-points are not enough to produce an accurate value for the MCA, we present in figure |2| the 
convergence of the MCA with the number of k-points for both compounds within LSDA. The MCA is the difference in 
the total energy between the in-plane axis [100] and [110] and the easy axis [001]. We have found that the total energy 
difference between the [100] and the [110] directions is negligible compared to the difference between the in-plane axis 
and the [001] axis. Our values are converged up to 6750 k-points and are given per unit cell (one atom of X and one 
of platinum). The MCA converges to 2.2 eV for CoPt and 3.9 eV for FePt. These behavior confirms the assumption 
that the system is isotropic inside the plane. The GGA MCA calculations converged to 1.9 meV and 4.1 meV for 
CoPt and FePt, respectively. The GGA results seem to be in good agreement with the LSDA results. 

Daalderop et a/o performed calculationsiJpr MCA in CoPt and FePtJpy means of an LMTO method in the atomic 
sphere approximation (ASA) within LSDAE3 using the force theorem.ell The easy magnetization axis found by this 
latter calculation is the [001] for both systems in agreement with our LSDA and GGA results. Their MCA value is 2 
meV for CoPt and is 3.5 meV for FePt. Solovyev et alxA used the same structure within a real-space Green's function 
technique framework to find also that the magnetization is along the [001] axis. Including the spin-orbit interaction 
for all the atoms, they found a CoPt MCA value of 2.3 meV and a FePt value of 3.4 ijneV. Our computed value 
agrees also with the value of 1.5 meV for CoPt and 2.8 meV for FePt, obtained by, Sakumac3 using LMTO method in 
the atomic sphere approximation (ASA) in conjunction with the force theorem.cil The drawback of the LMTO- ASA 
method is that it accounts only for the spherical part of the potential and ignores the interstitial region. Furthermore, 
the force theorem does not account directly for the exchange-correlation contribution to the MCA. Finally, Oppeneer 
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used the augmented spherical waves method (ASW) in the atomic sphere approximation and found a MCA value of 
2.8 meV for FePrtnand 1.0 meV for CoPt,c3 which are the smallest among all the ab-initio calculations. 

Grange et a/.Jlj using a torque measurement for a MBE deposited CoPt film on a MgO(OOl) substrate, obtained 
a MAE <?£il-0 meV. An early measurement of a monocristal of CoPt, by Eurin and Pauleve, produced a value of 
1.3 meVEa The large value of Eurin and_Paulcve is due to the fact that their sample was completely ordered. For 
FePt the first experiment of Ivanov et aZo produced an anisotropy value of 1.2 meV and showed that-for a thin film 
the^hapc anisotropy would be one order of magnitude smaller compared to the MCA. Farrow et alSB and Thiele et 
a/.Cj found for MBE deposited FePt films on a MgO(OOl) substrate an anisotropy value of 1.8 meV. These films were 
highly ordered (more than 95% of the atoms were in the correct site). All experiments have been carried out at room 
temperature, which explain at sope extend the difference between the calculated and experimental MCA values (the 
MCA decreases with temperatureEa). It is worth mentioning that the experimental MCA for FePt is much larger than 
that of CoPt in agreement with our calculations. For thick films volume shape anisotropy (VSA) contributes also to 
the MAE and it favors always an in-plane masaetization axis. We can estimate the VSA as — 27rMy- in c.g.s. units, 
where My is the mean magnetization densityj23 and obtain a value of -0.1 meV for FcPt and -0.06 meV for CoPt. 
These-values are one order of magnitude smaller than the MCA values in agreement with the speculation of Ivanov 
et alM 



IV. MAGNETIC MOMENTS 



A. Density of States 



In figure we present the cobalt-projected partial density of states for three spin quantization axis corresponding 
to angles 7= 0°, 45° and 90° calculated with 6750 k-points. The 3-c? states dominate the electronic structure of 
cobalt. The spin-up band is practically totally occupied while the spin down band is almost half-occupied. The 
general form of the cobalt DOS, as well as the iron DOS, does not seem to change appreciably with the angle 7. 
The platinum projected density of states show similar behavior for both FePt and CoPt compounds. Bulk platinum 
is paramagnetic and the small changes in the DOS come from the polarization of the 5d electrons via hybridization 
with the 3d electrons of cobalt (iron). It is worth mentioning that the DOS is calculated inside each muffin-tin and 
the interstitial region is not taken into account. To minimize the contribution of the interstitial region we use almost 
touching muffin-tin spheres. In addition, we have find that remaining region has a negligible spin polarization. 



B. Spin Magnetic Moments 



The spin magnetic moments are isotropic with respect to the spin quantization axis as expected. They were 
calculated by attributing all the charge inside each muffin-tin sphere to the atom located in that sphere. As outlined 
above we have found that the interstitial contribution to the spin magnetic moment is one order of magnitude smaller 
than that of the platinum site. In CoPt we found a cobalt spin magnetic moment of 1.79 and a platinum moment 
of 0.36 fiB within the LSDA. The GGA cobalt spin magnetic moment of 1.83 /i^ is slightly larger than the value 
within LSDA. This is due to a more atomic like description of the atoms in a solid within the GGA compared to 
the LSDA. Although the GGA underestimates the hybridization between cobalt and platinum d valence electrons 
compared to the LSDA, the larger cobalt spin moment leads to a slightly larger GGA platinum spin moment of 0.37 
fiB- Our calculated values are in good agreement with the experimental values of Grange et al. (1.75/iB for cobalt and 
0.35/iB for platinum).E3 Previous expetiments by van Laar on a powder sample gave a value of 1.7 /is for the cobalt 
atom and 0.25 for the platinum atom.c£l The spin magnetic mopients have been previously calculated by Solovyev et 
alxl (1.72 for cobalt and 0.37 for platinum), bjjJDaalderop et aO (1.86 for cobalt), by Sakumac3 (1.91 for cobalt and 
0.38 for platinum), and finally by Kootte et a/.cZl by means of a localized spherical wave method (1.69 for cobalt and 
0.37 for platinum). All previous calculations are in good agreement with our full-potential results. 

As expected the iron spin moments are much larger than the cobalt ones. The LSDA produced a value of 2.87 
fiB while the GGA produced a slightly larger value, 2.96 ^b- The hybridization between the iron 3d states and the 
platinum 5d states is less intense than in the case of CoPt resulting in a smaller platinum moment for FePt. The 
LSDA platinum spin moment is 0.33 ^b (compared to 0.36 /is in CoPt) and the GGA platinum spin moment is 0.34, 
/is (compared to 0.37 fiB in CoPt). The spin magnetic m|©inents have been previously caJculated by Solovyev et alcZI 
(2.77 for iron and 0.35 for platinum), by—Daalderop et aleJ (2.91 for iron), by SakumaeS (2.93 for iron and 0.33 for 
platinum), and finally by Osterloch et al£5 (2.92 for iron and 0.38 for platinum). Here again all previous calculations 
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are in good agreement with our full-potential results. All the methods produced a smaller induced spin moment for 
platinum atom in FePt than for CoPt, verifying that the hybridization effect in CoPt is much stronger than in FePt. 

C. Orbital Magnetic Moments 

Contrary to spin moments, the orbital moments are anisotropic. Figure ^ presents the behavior of the orbital 
moments as a function of the angle 7 between the [001] direction and the spin quantization axis within LSD A for both 
CoPt and FePt compounds (the lines are guide to the eye) . The orbital moments decrease with respect to the angle 
7 but the values for the [100] axis does not follow this general trend. All four lines seem to have the same behavior, 
but it is interesting to notice that for the CoPt and magnetization along the [100] axis cobalt orbital moment is 
smaller than the platinum one (for the values see Tables | and |^ . The cobalt moments decrease faster than platinum 
moments in CoPt. For the iron and platinum atoms in FePt the two lines are practically parallel. The platinum 
orbital moments in FePt are smaller than for the CoPt with a factor that varies from 73% for 7=0° to 67% for 7=90°. 
The ratio of iron and platinum orbital moments in FePt varies from 1.56 for 7=0° down to 1.35 for for 7=90°, which 
is considerably smaller than the ratio of cobalt and platinum orbital moments in CoPt. We notice here that in our 
calculations we can estimate only the projection of the total orbital moment on the spin quantization axis and we 
have no information concerning the real value of the total magnetic moment or its direction in space. It seems that for 
the magnetization along the [100] direction, the direction of the orbital magnetic moment undergoes a discontinuous 
jump, resulting in a large projection on the spin quantization axis with respect to 7 at the vicinity of 90°. This 
behavior is also reproduced by the GGA. 

In Table || we present the values of the orbital moments of iron and cobalt within both the LSDA and GGA as a 
function of the angle 7. The GGA values seem to be smaller than the LSDA ones but follow exactly the same trends. 
The orbital moment anisotropy is more important in the case of cobalt. The LSDA cobalt orbital moment changes 
by 0.048 fiB and the GGA moment by 0.027 fiB as we pass from the easy axis [001] to the hard axis [100]. The LSDA 
iron moment changes by 0.002 fiB and the GGA moment is the same for the two high symmetry directions. In the 
case of cobalt the difference between values calculated within the two functionals, LSDA and GGA, becomes smaller 
when the angle increases and for 75° it changes sign. In the case of iron this difference decreases only slightly with 
the angle but since the difference is considerably smaller than in CoPt (less than 0.004), we conclude that the orbital 
moment in the LSDA and the GGA are roughly the same. 

In Table ^ we present the values for the platinum orbital moment within both functionals. We see that the GGA 
produces larger moments than the LSDA for platinum in FePt contrary to CoPt. The platinum moments are in 
general smaller than the moments of the 3d ferromagnets, and the difference between the values calculated within 
LSDA and GGA are small. The absolute values for platinum are comparable to cobalt(iron) orbital moments even 
though the spin moments on platinum are one order of magnitude smaller than for cobalt (iron). The large orbital 
moments for platinum are due to a much larger spin-orbit coupling for the d electrons of the platinum compared to 
the 3d ferromagnets. |— ■ 

The orbital moments of £ePt and CoPt have been previously calculated by Daalderop and collaborators^ and by 
Solovyev and coUaboratorstZl for the [001] direction using the LSDA. The orbital moment of the cobalt site was found 
to be 0.12 iiB by Daalderop and 0.09 /is by Solovyev. The value of Daalderop is closer to our LSDA value of 0.11 
fXB- For iron site Daalderop found a value of 0.08 ^b and Solovyev 0.07 ^b in good agreement with our LSDA value. 
The platinum orbital moment has been calculated by Solovyev. He found a value of 0.06 ^b for platinum in CoPt 
and 0.044 /is for platinum in FePt, close to our values, 0.06 hb and 0.05 ^a- respectively. 

On the other hand, experimental data are available for CoPt by Grange;E3 obtained by applying the sum rules to 
the experimental XMCD spectra. The sum rules give the moments per hole in the d-band. To compare experiment 
with theory we calculated the number of d-holes by integrating the d projected density of states inside each muffin-tin 
sphere. We found 2.63 d and 2.48 d holes for cobalt and platinum, respectively. The cobalt orbital moment varies from 
0.26 ^B for 7=10° down to 0.11 /x_b for 7=60°. The measured values are also available for two other angles: 0.24 ^b 
for 7=30° and 0.17 fiB for 7=45°. Our theory reproduces qualitatively the experimental trends but underestimates 
the absolute values by more than 50% (see Tables | and |l| for all the values). The calculated values show a less sharp 
decrease with the angle than the experimental ones. For platinum the experimental data are available only for two 
angles 10° and 60°. For 7=10° the orbital moment is 0.09 /is and for 7=60° it is 0.06 /is. For the platinum site 
the calculated values are in much better agreement than for the cobalt site but we must keep in mind that the sum 
rules have their origin in an atomic theory and their use for 5d itinerant electrons like in platinum is still debated. 
The discrepancy between the theory and the experiment comes mainly from the approximation to the exchange and 
correlation. Both the LSDA and the GGA approximations to the density functional theory are known to underestimate 
the orbital moment values, because the orbital moment is a property directly associated with the current in the solid 
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and a static imag|e-|is not sufficient. But until now a DFT formalism like the current and spin density functional 
theory (CS-DFT)Jij which can treat at the same footing the Kohn-Sham and the Maxwell equations is too heavy to 
implement in a full-potential ab-initio method. The other problem is that there is no form of the exchange energy of a 
homogeneous electronic gas in a magnetic field known and this is the main quantity entering the CS-DFT formalism. 
Brooks has also developed an ad hoc correction to the Hamiltonian to account for the orbital polarization but this 
correction originates from an atomic theory and its application to itinerant systems is not satisfactory.tEl 

V. XMCD 

XMCD spectroscopy became popular after the development of the sum Pshlpj that enable the extraction of reliable 
information on the micro-magnetism directly from the experimental spectra.llSlla The great advantage of XMCD is that 
we can probe each atom and orbital in the system so to obtain information on the local magnetic properties. Lately 
angle-dependent XMCD experiments allowed the determination of magnetic ppaperties for different spin quantization 
axis. All experimental spectra for CoPt have been obtained by Grange et alji3 and all calculated spectra presented 
in this section were obtained using the LSDA. The GGA produced the same results and are not presented. 

Figure || presents the XMCD spectra for cobalt and iron atoms for the [100] and [001] magnetization axis. We 
convoluted_our theoretical spectra using a Lorentzian width of 0.9 eV and a Gaussian width of 0.4 eV as proposed 
by Ebert,E3 in the case of iron to account for the core hole effect and the experimental resolution, respectively. The 
energy difference between the L3 and L2 peaks is given by the spin-orbit splitting of the pi and pa core states. It 
is larger in the case of cobalt, 14.8 eV, than for iron, 12.5 eV. The intensities of the peaks are comparable for both 
atoms, but the cobalt peak-intensities are larger for the [001] axis contrary to the iron spectrum. The most important 
feature of these spectra is the integrated L3/L2 branching ratio as it enters the sum rules. It is larger in the case of 
cobalt site. This is expected since the sum rules predict that a larger L3/L2 ratio is equivalent to a larger orbital 
moment. For both atoms the integrated L3/L2 branching ratio is larger for the [001] axis. But the ratio anisotropy 
is larger for the cobalt site (1.32 for M||[001] and 1.17 for M||[100]) than for the iron site (1.15 for M||[001] and 1.14 
for M||[100]) reflecting the larger orbital moment anisotropy of cobalt compared to iron (see Table |). Especially 
for iron both orbital moment and integrated L3/L2 branching ratio are practically the same for both magnetization 
axis. As these changes in the XMCD signal are related to the change in the orbital magnetic moment between the 
two magnetization directions. The XMCD anisotropy should be roughly proportional to the underlying MCA but 
no relation exist thai connects these two anisotropies. However they are connected indirectly through the orbital 
moment anisotropyEa 

In figure^ we have plotted the platinum XMCD spectra for two angles 7=10° and 60° for both compounds. The life 
time of the core-hole in platinum is smaller than for cobalt so the broadening used to account for its life time should 
be larger. We used both a Lorentzian (1 eV) and a Gaussian (1 eV) to represent this life time and a Gaussian of 1 eV 
width for the experimental resolution. As in the case of cobalt and iron the platinum XMCD spectra change with the 
angle and depend on the surrounding neighbors. The peak intensities are larger in CoPt. Also the difference in the 
intensity due to the anisotropy has different sign in the two compounds. The intensities of FePt spectra for 7=10° 
are much larger than for 7=60° contrary to the CoPt behavior. As is the case for cobalt in CoPt, the platinum site 
integrated L3/L2 branching ratio shows larger anisotropy than for platinum in FePt. In CoPt it is 1.49 for 7=10° and 
1.19 for 7=60°, while for FePt it is 1.20 for 7=10° and 1.14 for 7=60°. Here again the XMCD follows the anisotropy 
of the orbital magnetic moment in these compounds. The energy difference between the two peaks is 1727 eV for 
both compounds. The 2p electrons of platinum are deep in energy and are little influenced by the local environment, 
so that their spin-orbit splitting does not depend on the neighboring atoms of platinum. 

We expect a better agreement between the theoretical and the experimental XMCD spectra for the platinum site 
than for the cobalt site, because the core hole is deeper and would effect less the final states of the photo-excited 
electron. In figure 0, we have plotted the absorption and the XMCD spectra of cobalt for 7=0°. We have scaled our 
spectra in a way that the experimental and theoretical L3 peaks in the absorption spectra have the same intensity. 
The energy difference between the L3 and L2 peaks is in good agreement with experiment. But the intensity of the 
L2 peak is larger than the corresponding experimental peak. The high intensity of the calculated L2 edge makes 
the theoretical XMCD integrated L3/L2 branching ratio of 1.32 much smaller than the experimental ratio of 1.72. 
This is because the LSDA fails to represent the physics of the core hole photo-excited electron recombination. In the 
case of 3d ferromagnets the core hole is shallow and influences the final states seen by the photo-excited electron. A 
formalism that can treat this electron-hole interaction has been proposed by Schwitalla and E^rt,E3 but it failed to 
improve the L3/L2 branching ratio of XAS of the late transition metals. Benedict and ShirleyE^I have also developed 
a scheme to treat this phenomena but its application is limited only to crystalline insulators. 

In figure ^ we have plotted experimental and theoretical total absorptions for the platinum atom in CoPt. For 
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both L2 and L3 edges the theory gives a sharp peak which does not exist in experiment. As expected the L3 peak is 
much more intense than the ^2- In contrast to what is obtained for cobalt, the results for the platinum XMCD (see 
figure ^) show better agreement with experiment, due to the fact that the core hole effect is less intense (core hole 
much deeper compared to cobalt). The experimental and theoretical L2 and L3 edges are separated by a spin-orbit 
splitting of the 2p core states of 1709 and 1727eV respectively. The width of both L2 and L3 edges is comparable to 
experiment, but the calculated L2 edge is much larger. This produces a calculated integrated branching ratio of 1.49 
which is much smaller than the experimental ratio of 2.66. Here again the theory is underestimating the branching 
ratio. 



VI. BAND STRUCTURE AND DENSITY OF STATES ANISOTROPY 

In figure |l^ we present the band structure along the [001] and [100] axis in the reciprocal space for different angles 
7 for the CoPt compound within the LSDA. We know that it is essentially the area around the Fermi level that 
changes with respect to the spin quantization axis. For this reason we have enlarged a region of ±1 eV around the 
Fermi level. In the first panel we plot the relativistic band structure for 7=0°. In the second and third panel we have 
plotted the relativistic band structure for 7=45° and 7=90°. We remark that as the angle increases there are bands 
that approach the Fermi level and cross it. However this information concerns just two high symmetry directions. 
For this reason we limited ourselves to the changes in the total DOS. In figure |ll| we notice that just below the Fermi 
level the DOS for the hard axis is lower than for the easy axis which seems to favor the hard axis. This means that 
the anisotropy does not originate from the changes ate the vicinity of the Fermi level but from that of the whole DOS. 
It is difficult to investigate this phenomena by inspection of the changes at the vicinity of t]a^ Fermi surface and to 
explain the sign of the MCA . Our results confirm the work of Daalderop and coUaboratorsEaJjhat argued that not 
only states in the vicinity of the Fermi surface contribute to the MAE, as originally thought j£3 but states far away 
make an equally important contribution. 

VII. CONCLUSION 

We have performed a theoretical ab-initio study of the magnetic properties of the ordered CoPt and FePt fct alloys 
systems. The calculated easy axis is the [001] for both compounds in agreement with other calculations and with 
experiments on films which found a strong perpendicular magnetization axis. The density of states is found to change 
very little with the direction of the spin quantization axis, and hence the magnetic moments are isotropic with respect 
to the magnetization axis. Contrary to the spin moments, the orbital magnetic moments decrease with the angle 7 
up to 75° (7 is the angle between the spin quantization axis and the [001] axis). 

The calculated x-ray magnetic circular dichroism (XMCD) for all the atoms reflect the behavior of the orbital 
moments. Especially platinum resolved spectra present large differences between the two compounds. Cobalt XMCD 
spectra are in agreement with experiment but as usual the L3/L2 ratio is underestimated by the theory. The platinum 
site shows better agreement with experiments because the core- hole is much deeper than in the case of cobalt. 

Finally we showed that all the occupied electronic states contribute to the magneto-crystalline anisotropy and not 
just states near the Fermi level for both CoPt and FePt. 
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TABLE L Calculated LSDA and GGA cobalt(iron) orbital magnetic moments with respect to the angle 7 between the [001] 
axis and the spin quantization axis in the (010) plane. For both atoms the orbital moments decrease with the angle but values 
for 7=90° do not follow this trend. Both LSDA and GGA produce the same trends. The moments at the cobalt site are 
larger than at the iron site. The experimental orbital moments for the cobalt site are taken from the Ref. |l9| The theory 
underestimates the experimental orbital moments by more than 50% but reproduces the correct trends. 
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TABLE II. Calculated LSDA and GGA orbital moments of platinum in FePt and CoPt with respect to the angle 7 between 
the [001] axis and the spin quantization axis in the (010) plane. For the FePt the GGA produces larger values than LSDA 
contrary to CoPt. The moments decrease with the angle but increase sharply at 7=90°. The moments are significantly larger 
for the CoPt compound. Experimental results for platinum in CoPt are taken from the Ref. |l^. The agreement is better than 
in the case of cobalt. 
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FIG. 1. Total energy for the CoPt system within LSDA as a function of the angle 7 between the [001] axis and the spin 
quantization axis on the (010) plane. The fundamental state corresponds to 7=0°, that means to the [001] axis. The total 
energy increases as a function of 7. 
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FIG. 2. Convergence of the magneto crystalline anisotropy in meV with respect to the number of k-points. The calculation 
is performed for two in-plane axis. For CoPt the line for the [100] axis converges to 2.2 meV while the line for the [110] axis 
converges to 2.17 mcV. For FePt the two lines differ by less than 0.01 meV. The behavior of MCA verifies the assumption that 
the Llo structure is isotropic in the (001) plane. 



10 




-10 -5 

Energy (eV) 



FIG. 3. Cobalt projected partial density of states with respect to the angle 7 between the [001] axis and the spin quantization 
axis in the (010) plane. There are no major changes with 7. This behavior reflect the isotropic character of the spin moments. 
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FIG. 4. Calculated orbital moments for both compounds within LSDA as a function of the angle 7 between the [001] 
direction and the spin quantization axis (the lines are guides for the eye). All four atoms show the same behavior. The orbital 
moments decrease with the angle 7 until about 75°, then they increase. The moments for the CoPt compound are larger than 
for these of FePt. 
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FIG. 5. XMCD spectra for the cobalt and iron sites in the two XPt compounds and for magnetization along the two high 
symmetry axis ([001] and [100]). The iron and cobalt XMCD show different behavior. The picks intensities for the cobalt 

spectrum lay higher for M||[001] and the L-^/L2 integrated branching ratio is larger for this axis. Contrary, the picks intensity 
for the iron site are larger for M||[100] and the L3/L2 integrated branching ratio is about constant. 
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FIG. 6. XMCD spectra for platinum for two values of the angle 7 between the [001] direction and the spin quantization 
axis in both compounds. The differences in the XMCD are important and the L3/L2 integrated branching ratio changes 

considerably with the angle. The intensities for the CoPt compound arc larger reflecting the bigger orbital moments in this 
compound. The spin-orbit splitting of the core pi and pa states is 1727eV in both compounds. 
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FIG. 7. Theoretical and experimental absorption and XMCD spectra for the cobalt site in CoPt for 7=0°, where 7 the 
angle between the [001] direction and the spin quantization axis. The theory overestimates the absorption at the L2 edge and 
underestimates the L3 / L2 integrated branching ratio. 
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FIG. 8. Calculated and experimental absorption spectra at the platinum site and for 7=10° in CoPt, where 7 is the angle 
between the [001] direction and the spin quantization axis. The theory produces a peak at the threshold of both L2,3 edges. 
Both experimental and theoretical spectra have the same shape for the two edges but L3 edge is much more intense. Most of 
the structures in the experimental spectra are washed out due to a strong broadening effect. 
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FIG. 9. Calculated and experimental platinum XMCD spectra at L2 and L3 edges for 7=10° in the CoPt compound (7 is 
the angle between the [001] direction and the spin quantization axis). The agreement between the theory and the experiment 
is better than in the case of cobalt because the core hole is deeper. 
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FIG. 10. CoPt band structure around the Fermi level along two high symmetry directions of the Brillouin zone as a function 
of the angle 7 between the [001] direction and the spin quantization axis. 
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FIG. 11. Total density of states for CoPt at the vicinity of the Fermi level, for the spin quantization axis along the [001] and 
[100] magnetization axis. The states just below the Fermi level suggest that the hard axis is favored which led as to conclude 
that the MCA is due to the changes in all the DOS. 
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